Mon. Not. R. Astron. Soc. 000, ITHTT1 (20131 Printed 11 March 2013 



(MN WFeX. style file v2.2) 



A near-infrared relationship for estimating black hole 
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ABSTRACT 

Black hole masses for samples of active galactic nuclei (AGN) are currently esti- 
mated from single-epoch optical spectra using scaling relations anchored in reverber- 
ation mapping results. In particular, the two quantities needed for calculating black 
hole masses, namely, the velocity and the radial distance of the orbiting gas are de- 
rived from the widths of the Balmer hydrogen broad emission lines and the optical 
continuum luminosity, respectively. We have recently presented a near-infrared (near- 
IR) relationship for estimating AGN black hole masses based on the widths of the 
Paschen hydrogen broad emission lines and the total 1 /iin continuum luminosity. The 
near-IR offers several advantages over the optical: it suffers less from dust extinction, 
the AGN continuum is observed only weakly contaminated by the host galaxy and the 
strongest Paschen broad emission lines Paa and Pa/3 are unblended. Here we improve 
the calibration of the near-IR black hole mass relationship by increasing the sample 
from 14 to 23 reverberation-mapped AGN using additional spectroscopy obtained with 
the Gemini Near- Infrared Spectrograph (GNIRS). The additional sample improves the 
number statistics in particular at the high luminosity end. 
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1 INTRODUCTION 

The discovery of tight correlations between a galaxy's 
central black hole mass and the luminosit y, velocity dis- 
persi o n and mass of i ts ste ll ar bulge jMagorrian et al.l 



19981: iGebhardt et all |2000|; [Ferrarese fc Merrittl l200d: 



Marconi fe Hunrj|200j : lHaring fc~R ix 2004) is expected to 



strongly constrain how galaxies form and grow over cosmic 
time. Therefore, much effort in particular goes into mea- 
suring the rate of black hole growth (e.g.. | Yu fc Tremainei 
120021 : iHeckman et alJliool : iKellv et al.ll2009l ). Since this re- 
quires both large samples of galaxies with easily obtainable 
black hole mass estimates and sources that probe the highest 
redshifts, such studies rely heavily on active galactic nuclei 
(AGN). 

In general, the masses of black holes at the centres of ac- 
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tive galaxies are measured based on the gravitational force 
they exert on other massive objects such as stars or gas 
clouds. But, whereas measurements of the stellar dispersion 
require spectroscopy at very high angular resolution in or- 
der to resolve the black hole's sphere of influence, the ve- 
locity dispersion of the broad-emission line gas present at 
the centre of AGN can be readily measured from a single 
long-slit optical, ultraviolet (UV) or near-infrared (near-IR) 
spectrum. Then, assuming that the dynamics of the broad- 
emission line gas is dominated by gravitational forces, one 
can use the virial theorem to calculate black hole masses: 



M BH = / 



RAV 2 



(1) 



where R is the radial distance of the broad emission line gas 
from the black hole, AV is the velocity dispersion of the gas, 
G is the gravitational constant and / is a scaling factor that 
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depends on the (unknown) dynamics and geometry of the 
broad line region (BLR). 

The BLR radius can be directly measured through 
reverberation mapping, a technique which determines the 
light-travel time-delayed lag with which the flux of the 
BLR responds t o changes in th e ionising continuum flux 
(see reviews by iPetersonl 1 19931 ; iNetzer fc Peterson! 1 19971 ; 
iPeterson fc Bent j|2006h ~ However, since reverberation map- 
ping campaigns are observing time intensive, at present 
the BLR radii of fewer than 50 A GN have been measured 
i Kaspi ct al. 2000; IPeterson et al.l 



iBentz et al.l |200S| ; iBentz et al.1 l2009bl ; iBarth et al.l l201lf ) 



2004; iKaspi et al 



2007; 



Therefore, the so-called radius-luminosity (R-L) relation- 
ship is u sed to est i mate t he BLR sizes for lar g e samples of 
AGN As |Petersod(ll993l).IWandel et al.l (| 19991 ) . IKaspi et al.1 
( 2000) and IBentz et al.1 i|2009al ) have shown, the BLR lags 
obtained from reverberation mapping campaigns correlate 
with the optical luminosity (of the ionising component) 
largely as expected from simple photoionisation arguments. 
The continuum luminosity L is usually measured in the op- 
tical (at rest-frame 5100 A) or in the UV. 

Given the importance of the R-L relationship for AGN 
black hole mass determinations, alternatives to the optical 
and UV continuum luminosity are being explo red, e.g., the 
X-ray luminosity and broa d H/3 line luminosity (|Kaspi et al.l 
l2005l ; iGreene et all l201fj|). the JOiii] A5007 A emission- 



line luminosity ( Greene et al.l |2010|) , the Paa and Pa/3 
emission-line luminosities ( Kim et al] |2010). and the mid- 
infrared [Ne v] A 14.32 ^m and [Q iv] A25.89 nm emiss ion- 
line luminos it ies foasyra et al . 2008; G reene et alj|201oh . In 
lLandt et al.l (|2011al ) we proposed the 1 /jm continuum lu- 
minosity as an efficient alternative that connec ts directly to 
the spe ctrum of the ionising source. Then, in lLandt et al.l 
|2011bl ) we introduced the relationship between the black 
hole mass and the near-IR virial product, i.e., the product 
between the total 1 /jm continuum luminosity and the width 
of the Paschen hydrogen broad emission lines. 

The near-IR potentially affords some advantages over 
the optical and UV. First, the observed optical continuum 
can be severely contaminated by host galaxy starlight if a 
large slit is used (as is often the case in reverberation stud- 
ies), especially in low-redshift sources that have a weak AGN 
or sources with a luminous stellar bulge, and so does not give 
directly the ionising flux. High-resolution, deep optical im- 
ages are then required to estimate the host galaxy starlight 
encl osed in the aperture in orde r to correct the optical spec- 
tra jBentz et al.l l2006bl . l2009al ). Secondly, the optical hy- 
drogen broad-emission lines, in particular H/3, and most of 
the UV broad emission lines are strongly blended with other 
species, making it difficult to reliably measure the line width. 
Thirdly, all optical and UV measurements may suffer from 
dust extinction. 

Near-IR spectroscopy is now available at excellent ob- 
serving sites that regularly achieve subarcsecond seeing. Ad- 
ditionally, the seeing is improved at longer wavelengths (e.g., 
by a factor of ~ 1.3 at 1 /im relative to 5100 A). This means 
that the host galaxy flux contribution in the near-IR can 
be reduced to a minimum. Note, however, that for the same 
(small) aperture the optical spectrum would be less domi- 
nated by host galaxy light than the near-IR one since the 
AGN to host galaxy flux ratio is intrinsically higher at op- 
tical wavelengths. Then, as we have shown in lLandt et al.l 
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Figure 1. The radius of the H/3 broad-emission line region (in 
light-days) versus the black hole mass (in solar masses) for the 
reverberation-mapped AGN sample with (filled circles) and with- 
out near-IR spectroscopy (open circles). 



(2008), the strongest Paschen hydrogen broad emission lines 
Paa and Pa/3 are observed to be free from strong blends 
offering the opportunity to reliably measure their widths. 
Furthermore, since attenuation by dust is much reduced in 
the near-IR compared to the optical and UV, our near-IR 
relationship should be particularly useful as an alternative 
for dust-obscured AGN. 

Here we present an improved calibration of our near- 
IR relationship for estimating black hole masses in AGN. 
The paper is structured as follows. In Section 2, we intro- 
duce the enlarged sample that includes additional sources, 
which were observed as described in Section 3. In Section 4, 
we present an improved calibration of the near-IR R-L rela- 
tionship, and in Section 5, we derive the revised relationship 
between the black hole mass and the near-IR virial product. 
A few sources in our sample are found to have near-IR lumi- 
nosities dominated by host galaxy starlight and we discuss 
them in the context of the near-IR relationships in Section 
6. In Section 7, we summarize our results and present our 
conclusions. Throughout this paper we have assumed cos- 
mological parameters Ho = 70 km s _1 Mpc -1 , Dm = 0.3, 
and Qa = 0.7. 



2 THE SAMPLE 



IPeterson et al.l (|2004T l reanalysed the available reverberation 
mapping data of a sample of 37 AGN and presented im- 
proved black hole m as s determination s for 3 5 objects. Later, 
IBentz et al.1 ll2006al). Ibenney et al.1 (|2006l ), iDennev et all 
( 20101 ) and iGrier et al.l (|2012l ) refined the black hole mass 
estimates for some of these sources and added a new deter- 
mination for the source Mrk 290. 

Ou r original sample o bserved with the SpeX spectro- 
graph (|Ravner et al.l 120031 ) on the NASA Infrared Tele- 
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Table 1 . General properties of the additional sample 



Object name 
(1) 


R.A. (J2000.0) 

(2) 


Decl. (J2000.0) 
(3) 


z 

(4) 


J 

(mag) 
(5) 


Rn/3 
(lt-days) 
(6) 


Ref. 

(7) 


A/bh 
(10 6 M ) 
(8) 


Ref. 
(9) 


log vLl^m 

(erg s _1 ) 
(10) 


line 
(11) 


type 
(12) 


FWHM 
(km s" 1 ) 
(13) 


""line 

(km s _1 ) 
(14) 


PG 0026+129 


00 29 13.60 


+13 16 03.0 


0.145* 


13.58 


"US 


B09 


393+96 


P04 


44.72 


Pact 


e 


1793 


1798 


PG 0052+251 


00 54 52.10 


+25 25 38.0 


0.155 


13.89 


B09 


369+76 


P04 


44.46 


Pact 


i 


4133 


1750 


3C 120 


04 33 11.10 


+05 21 15.6 


0.033 


12.67 




G12 


66+7 


G12 


43.45 


Pa^ 


i 


2733 


1815 


PG 0804+761 


08 10 58.60 


+76 02 42.5 


0.100 


12.98 


147+19 


B09 


693+83 


P04 


44.58 


Pact 


lack 


2304 


1193 


NGC 3516 


11 06 47.49 


+72 34 06.9 


0.009 


11.13 


12±1 


D10 


32±l 


D10 


43.04 


Pa^ 


lack 


4469 


1527 


NGC 4051 


12 03 09.61 


+44 31 52.8 


0.002 


11.65 


1.9+0.5 


D10 


1.7+°'° 

— . 5 


D10 


41.51 


Pa,3 


e 


1681 


762 


PG 1211+143 


12 14 17.70 


+14 03 12.6 


0.081 


13.15 


Q4+26 


B09 


146+44 


P04 


44.33 


Pact 


c 


1601 


955 


PG 1307+085 


13 09 47.00 


+08 19 48.2 


0.155 


14.16 


106+^6 


B09 


440+123 


P04 


44.61 


Pact 


i 


2982 


1631 


Mrk 279 


13 53 03.45 


+69 18 29.6 


0.030 


12.19 


17+4 


B09 


35+9 


P04 


43.43 


Pa^ 


i 


3568 


1580 


3C 390.3 


18 42 09.00 


+79 46 17.1 


0.056 


13.60 


24+« 


B09 


287+64 


P04 


44.03 


Pact 


i 


7460 


2612 


Mrk 1513 


21 32 27.81 


+10 08 19.5 


0.063 


12.84 


n n+0.6 

LO -°-0.5 


G12 


41+9 


G12 


44.08 


Pact 


e 


1905 


1265 



The columns are: (1) object name; ( 2) and (3) positio n, and (4) redshift from the NASA/IPAC Extragalactic Database (NED); (5) J band magnitude from the 
Two Micron AH- Sk y Survey C2MASS; ICutri et al 1120031); (61 radius of t he H/3 broad -emission line region (in light-days); (7) reference for the H/3 radius, where B09: 
Bcnt z et al. 1 1 1200931 ), D10: lDennev et al.l j201Ch . G12: lGrier et alj J2012T); (8) black hole mass (i n solar mass e s) calculated from the line dispers ion assuming a scaling 
factor of / = 5.5; (9) reference for the black hole mass, where DlOi lDennev et al.l l|2010h . G12: iGrier et alj j2012h . P04: IPeterson et al] \2004 ); (10) integrated total 
l-/im continuum luminosity; (11) Paschen emission line with the highest signal-to-noise ratio; (12) type of transition between the broad and narrow emission line 
component, where i: inflected (i.e., transition point is obsvious), e: estimated (i.e., transition point is estimated) and lack: absent narrow emission line component; 
(13) full width at half maximum (FWHM) and (14) dispersion of the broad emission line component (both uncorrected for instrumental resolution). Errors represent 
Iff uncertainties. 

* We measure a redshift of z = 0.145 from narrow emission lines, whereas a value of z = 0.142 is listed in NED. 
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NGC 3516 (z = 0.009) 
2 1 0.8 0.6 0.4 



NGC 4051 (z = 0.002) 
2 1 0.8 0.6 0.4 





PG 0052 + 251 (z = 0.155) 



Mrk 279 (z = 0.030) 



aO 
o 





l0 g ^rest [ HZ 1 

Figure 2. Rest-frame spectral energy distributions (SEDs) for the two sources (NGC 3516 and NGC 4051) that are strongly affected 
by host galaxy light (top panels) and for the two sources (PG 0052+251 and Mrk 279) that have an SED dominated by AGN emission 
although estimates based an HST image decomposition give relatively high host galaxy fl ux contributions (bott om panels). The near-IR 
spectr a are normalized at 1 /an (dashed line) and non-simultaneous optical spectra f rom iKe nnicuttl dl992Tl and iMoustakas fe Kennicuttl 
(2006) arc also shown for NGC 3516 and NGC 4051. The host galaxy templates from lMannucci et al. are overplotted (magenta). 

Wavelength units in (im are labeled on the top axis. 



scope Facility (IRTF), a 3 m telescope on Mauna Kea, 
Hawai'i, included 16/36 re verberation-mapp ed AGN from 
the Peterson et al. sample (|Landt et al.ll2008l ). For observa- 
tions with the Gemini Near-Infrared Spectrograph (GNIRS) 
at the Gemini North observatory, an 8.1 m telescope on 
Mauna Kea, Hawai'i, we have selected all the remaining 
northern (S > —37°) reverberation-mapped AGN. We were 
granted observations of one source (3C 120) in the Sci- 
ence Verification (SV) phase of GNIRS and of eight sources 
in semester 2011B. Additionally, two reverberation-mapped 
A GN (NGC 4051 and Mrk 279) were observed with the IRTF 
by iRiffel et al l (120061 ) and we made use of their publicly 
available spectra. Therefore, our enlarged sample contains 
27/36 reverberation-mapped AGN from the Peterson et al. 
sample. Table Q] lists the general properties of the additional 
sample. 

Fig. Q]shows the radius of the H/3 broad-emission line re- 
gion (Rhp) versus the black hole mass for the reverberation- 
mapped AGN with (filled circles) and without near-IR 
spectroscopy (open circles). The remaining sample of nine 
sources would increase the number statistics in particular at 
the high end of values. Therefore, obtaining near-IR spec- 



troscopy of them could lead to a significant improvement 
of the calibration of our near-IR relationship for estimating 
black hole masses. 

Four sources in our enlarged sample, namely, Mrk 590, 
NGC 3227, NGC 3516 and NGC 4051, were found to be in 
a very low luminosity state with a continuum spectral en- 
ergy distribution (SED) strongly dom inated by ho s t galax y 
emission (see Fig. top panels, and lLandt et all (|2011bh . 
their Fig. 6). We have excluded these sources from the cali- 
bration of the near-IR relationships in Sections [4] and [5] but 
we discuss them in the general context in Section^ since an 
upper limit to the integrated 1 /im continuum luminosity of 
the AGN can be obtained from the total observed value. 



3 THE OBSERVATIONS 

We observed in queue mode with the G emini Near-Infrared 
Spectrograph (GNIRS; lElias et al1l2006l ). which has recently 
been recommissioned at the Gemini North observatory. One 
source (3C 120) was observed during the Science Veri- 
fication (SV) phase of GNIRS in semester 2010B (Pro- 
gram ID: GN-2010B-SV-165) and eight sources in semester 
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Table 2. Gemini Journal of observations 



Object name 


observation 


exposure 


airmass 


continuum S/N 


telluric standard star 




date 


(sec) 




J 


H 


K 


name 


distance 


airmass 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(°) 
(9) 


(10) 


PG 0026+129 


2011 Aug 3 


8x120 


1.144 


62 


95 


69 


HD222558 


12.1 


1.113 


PG 0052+251 


2011 Aug 3 


8x120 


1.012 


69 


89 


119 


HD222558 


23.1 


1.149 


3C 120 


2010 Dec 15 


2x120 


1.046 


49 


55 


164 


HD34035 


10.7 


1.053 


PG 0804+761 


2011 Nov 30 


8x120 


1.828 


84 


106 


160 


HD48049 


8.8 


1.786 


NGC 3516 


2011 Dec 8 


4x120 


1.791 


130 


168 


149 


HD48049 


20.2 


1.606 


PG 1211+143 


2011 Dec 18 


8x120 


1.473 


147 


70 


119 


HD101060 


8.9 


1.402 


PG 1307+085 


2011 Aug 11 


8x120 


1.851 


50 


48 


92 


HD116960 


5.5 


1.789 


3C 390.3 


2011 Aug 4 


8x120 


2.034 


132 


156 


187 


HIP942 


27.2 


1.901 


Mrk 1513 


2011 Aug 3 


8x120 


1.034 


132 


92 


85 


HD210501 


11.9 


1.116 



The columns are: (1) object name; (2) observation date; (3) exposure time; (4) average airmass; S/N in the 
continuum over ~ 100 A measured at the central wavelength of the (5) J, (6) H, and (7) K band; for the star 
used to correct for telluric absorption (8) name, (9) distance from the source, and (10) average airmass. 



2011B (Program ID: GN-2011B-Q-97). We used the cross- 
dispersed mode with the short camera at the lowest res- 
olution (31.7 1 mm -1 grating), thus covering the entire 
wavelength range of 0.9 — 2.5 (im without inter-order con- 
tamination. We chose a slit of 0.3" x 7" for the source 
3C 120 and, due to a poor weather allocation, a larger slit 
of 0.675" x 7" for the sources observed in semester 2011B. 
This set-up gives an average spectral resolution of full with 
at half-maximum (FWHM) ~ 180 km s _1 for 3C 120 and 
of FWHM~ 400 km s~ for the remainder of the sources. 
The chosen exposure times ensured that we obtained spec- 
tra with a high signal-to- noise ratio (S/N ~ 100) in order to 
reliably measure the broad emission line profiles. 

Either before or after each source, we observed a nearby 
(in position and air mass) AO V star with accurate near- 
IR magnitudes. These stars were used to correct the source 
spectra for telluric absorption and for flux calibration. Flats 
and arcs were taken with each source/telluric standard star 
pair. In Table [5] we list the journal of observations. The 
data were reduced using the Gemi ni/IRAF package (versio n 
1.11) with GNIRS specific tools (|Cooke fc Rodgersll2005l ). 
The data reduction steps included preparation of calibra- 
tion and science frames, processing and extraction of spectra 
from science frames, wavelength calibration of spectra, tel- 
luric correction and flux-calibration of spectra, and merging 
of the different orders into a single, continuous spectrum. 
The spectral extraction width was adjusted interactively for 
each telluric standard star and source to include all the flux 
in the spectral trace. The final spectra were corrected for 
Galactic extinction using the IRAF task onedspec.deredden 
with input Av values deri ved from Galactic hydroge n col- 
umn densities published bv lDickev fc Lockmanl (|1990l ). The 
results are shown in Fig. IA1I 



4 THE NEAR-IR R-L RELATIONSHIP 

In lLandt et all (|2011bh we have shown that the accretion 
disc spectrum of AGN, which is believed to be the main 
source of ionising radiation, extends well into the near-IR 
and still dominates the continuum at ~ 1 /im. Therefore, 
a single-epoch near-IR spectrum of a broad-emission line 
AGN can in principle be used to estimate the BLR radius. 




Figure 3. The radius of the H/3 broad-emission line region (in 
light-days) versus the integrated total 1 /im continuum luminos- 
ity. Black and red filled circles indicate the original sample of 
lLandt et al. I ll2011at) and our additional sample, respectively. The 
sources with a continuum SED dominated by host galaxy emis- 
sion are shown as green filled circles. The dashed, solid and dotted 
lines indicate the correlations obtained with the BCES, FITEXY 
and GaussFit routines, respectively. 



We verified this conjecture in lLandt et al.l (|2011al ). where 
we introduced the near-IR radius - luminosity (R-L) rela- 
tionship. Here we present an improved calibration of it. 

In Fig. we plot the radius of the H/3 broad-emission 
line region (Rhp) versus the integrated total 1 fim contin- 
uum luminosity for our enlarged sample of reverberation- 
mapped AGN (27 objects). Value s for the r a dius o f 
the H/3 broa d -line regio n are from [Bentz et all (|2009al ) . 
iDennev et ail (|2010l ) and iGrier et al l (|2012l ) (see also Ta- 
ble [TJ. The near-IR measures for the original sample (black 
filled circles) are based on our IRTF data, whereas we have 
measured the integrated total 1 fim continuum luminosity 
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Table 3. Estimates of the near-IR host galaxy contribution 



Object Name 


a 

(1 + 2)5100 


host 


rcf . 


near-IR spectrum 


HST flux 


lop; v Tj-\ ,, m 


AGN/ 




( map") 


tvnp 




aperture 
(arcsec 2 ) 


PA 

(°) 


((1 + z)5100 A) 
(erg s ^ cm 2 A ^) 


(prp- c{ 1) 


host 


en 


(2) 


(3) 




(5) 


(6) 


(71 


(81 


(9) 




0.105 




B09 


0.7 x 5.8 





"7 a k T? 1 1 
1 .bobj—L / 


43.70 


9 


PG 0052+251 


0.115 


Sb 


H02 


0.7 x 4.3 





2.67E-16 


44.22 


0.7 


3C 120 


0.509 


SO 


P07 


0.3 x 1.9 





7.30E-17 


42.26 


14 


PG 0804+761 


0.014 


E 


B09 


0.7 x 5.5 





2.42E-17 


42.83 


55 


NGC 3516 


0.030 


Sb 


RC3 


0.7 x 4.5 





3.06E-15 


42.66 


1 


NGC 4051 


0.000 


Sb 


RC3 


0.8 x 1.6 


132 


2.04E-15 


41.17 


1 


PG 1211+143 


0.000 


E 


B09 


0.7 x 5.5 





1.63E-16 


43.46 


6 


PG 1307+085 


0.000 


E 


B09 


0.7 x 3.7 





2.94E-17 


43.35 


17 


Mrk 279 


0.000 


Sa 


P07 


0.8 x 1.7 





1.12E-15 


43.33 


0.3 


3C 390.3 


0.085 


Sa 


B09 


0.7 x 5.6 





3.58E-16 


43.40 


3 


Mrk 1513 


0.114 


SO/a 


P07 


0.7 x 5.7 





1.68E-16 


43.22 


6 



The columns are: (1) object name; (2) Galactic ext inction at rest-f r ame 5 100 A; (3) Hubble typ e of th e host 
galaxy; (4) referenc e for t he host type, where H02: lHamilton et al.l ll2002h . P07: iPetrosian et all (1200711 . B09: 
IBentz et al]||2009all , RC3: lde Vaucouleurs et~al"1 lll99ll) ; (51 near-IR extraction aperture; (6) near-IR slit position 
angle, where PA=0° corresponds to E-W orientation and is defined E through N; (7) host galaxy flux at rest- 
frame 5100 A, corrected for Galactic extinction using the values in column (3); (8) integrated host galaxy l-/im 
continuum luminosity; and (9) luminosity ratio between AGN and host galaxy at rest-frame l-fim. 



Table 4. Best-fits for the relation log_R H /3 = K + a \og{uL v /10 4z ) 





original sample (14 obj.) 


enlarged sample (23 obj.) 


Type 


K 


a 


S* 


K 


a 


S* 


BCES 


Avg 
MC 


0.91+0.11 
0.90+0.13 


0.56+0.10 

56+° 12 




0.89+0.11 
0.88+0.12 


0.64+0.09 
0.64+0.10 




FITEXY 


Avg 
MC 


0.82+0.10 
0.81+0.11 


0.52+0.10 
0.53+0.11 


53 
52+2 


0.83+0.09 
0.81+0.10 


0.57+0.08 
0.58+0.09 


53 
52+1 


GaussFit 


Avg 
MC 


1.00+0.10 
0.95+0.14 


0.54+0.10 
0.56+0.12 




0.94+0.09 
91+ 012 


0.60+0.09 
0.62+0.11 





* Scatter calculated as the percentage of the logijjj^ value that, when added 
in quadrature to the error value, gives = !• 



for the additional sample (red filled circles) in the GNIRS 
spectra obtained at the Gemin i North observator y and for 
Mrk 279 in the IRTF data of iRiffel et ail (|2006t ). We ob- 
served the original IR TF sample on ave r age tw ice within a 
period of ~ 3 years l|Landt et all I2008I . bOllbl ). For these 
sources we use the mean integrated total 1 /im continuum 
luminosity and the error on the mean. For the Gemini sam- 
ple, which has only one observation epoch available, we have 
adopted the average error of the IRTF sample of 0.04 dex. 

The excellent atmospheric seeing on Mauna Kea, 
Hawai'i, where both the IRTF and the Gemini North ob- 
servatories are located, allowed us to use a relatively narrow 
slit (< 1"), which excluded most of the host galaxy starlight. 
We have verified that our spectra are indeed dominated by 
the AGN in two separate ways. First, we have estimated 
the host galaxy contribution in the near-IR ap erture using 
the Hubble Space Telescope (HST) images of iBentz et all 



(|2006bh and IBentz et all (|2009al ) and following their ap- 
proach. The observed HST fluxes were transformed to a 
rest-frame wavelength of 5100 A by applying a colour cor- 
rection factor based on the model bulge galaxy template of 
iKinnev et ail l|l99fj) and were corrected for Galactic extinc- 
tion usin g A\ values derived f rom th e hydrogen column den- 
sities of Dickey fc Lockmanl (1990). The unabsorbed rest- 
frame 5100 A fluxes were then used to derive rest-frame 1 fim 
fluxes by scaling the galaxy template from iMannucci et all 
l|200lf) of the appropriate Hubble type. Details of our host 
galaxy flux estimat es for the original sample are listed in 
lLandt et all (|2011tJ l and for the additional sample are given 
in Table [3] Secondly, since HST image decomposition of 
AGN can be problematic at the centre of the galaxy where 
the bright AGN is degenerate with the concentration of the 
compact bulge, we have checked that the shape of the ob- 
served continuum SED is dominated by the AGN, i.e. that 
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we observe the accretion disc and the hot dust emission ris- 
ing and falling towards longer frequencies (in a logarithmic 
vL v versus v plot), respectively, resulting in a characteris- 
tic inflection point at the location where they meet (at ~ 1 
/an). 

As expected, the relatively small aperture used in the 
near-IR leads to relatively high AGN to host galaxy flux ra- 
tios and so we observe a continuum SED dominated by the 
AGN. However, four sources (NGC 35f6 and NGC 4051, 
see Fig. (2) top panels, an d Mrk 590 and NGC 3227, see Fig. 
6 in iLandt et all l|2011bl )) were found to have a continuum 
SED strongly dominated by host galaxy emission. Since in 
these cases the observed 1 /tm continuum luminosity rep- 
resents only an upper limit for the AGN contribution, we 
have excluded these from the following analysis. For two 
sources (PG 0052+251 and Mrk 279) the decomposition of 
HST images gives a relatively high host galaxy flux contri- 
bution, however, their continuum SED is clearly dominated 
by AGN emission (see Fig. [2l bottom panels) . 

With the additional sample (red filled circles) we have 
improved the number statistics in particular at the high 
end of the R—L relationship, where we previously had 
only two sources. More than half of the additional AGN 
have radii of the H/3 broad-emission line region of Rug > 
50 light-days and total 1 /im continuum luminosities of 
vLipm J> 2 ■ 10 44 erg s _1 . At the low-luminosity end we 
find mainly sources with a continuum dominated by host 
galaxy light (green filled circles). Therefore, improving the 
number statistics of the R-L relationship in this regime and 
extending its dynamic range to even lower values is likely to 
be problematic. 

We have performed linear fits of the form 
lo g Rub = K + a log ( vL v / 1 43 ) followin g the approach 
of iBentz et all (|2006bh and iBentz et all (|2009al ). In par- 
ticular, we have used the three fitti ng routines BCES 
l|Akritas fc Bershadvl 1 19961). FITEX Y i|Press et alj I2007T I 
and GaussFit i McArthur et al.lll994l ) that can incorporate 
errors on both variables and, except for GaussFit, allow us 
to account for intrinsic scatter. Note that accounting for 
intrinsic scatter has the effect of increasing the weight given 
to data points with the largest errors, which is preferred 
if the i ntrinsic dispersion is larger than the measurement 
errors (|Tremaine et all 120021 '). Roughly a quarter of our 
sample (6/23 sources) has multiple measurements of Rug 
and, therefore, we have considered the following two cases: 
(i) using the average derived from all measurements for a 
particular source and weighted by the mean of the positive 
and negative errors (i.e., weighted averages) and (ii) using 
Monte-Carlo (MC) techniques to randomly sample Rhb 
from the individual values for each o bject. In Ta ble [4] we 
compare our previous results from ILandt et al.l (|2011al ) 
with the results for the enlarged sample. The latter are 
shown for the case of the weighted averages in Fig. [3] Note 
that the near-IR continuum luminosities and the optical 
broad-emission line radii are not measured simultaneously, 
which may increase the scatter in their relationship. 

Simple photoionisation arguments suggest that a given 
emission line will be produced at the same ionising flux in all 
AGN and therefore R oc L 1//2 . We obtain with the enlarged 
sample in all six cases a best-fitting slope of ~ 0.6 ±0.1, 
which is slightly steeper than our previous results but con- 
sistent with a value of 0.5. The errors on the slope and 
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Figure 4. Black hole mass determined from optical reverberation 
campaigns versus the near-IR virial product between the total 
continuum luminosity at 1 /im and the line width of the Pact or 
Pa/3 broad component, whichever line had the higher signal-to- 
noise ratio (left panel) and versus the black hole mass calculated 
using the near-IR R-L relationship and the virial theorem (right 
panel). We have used for the line width the full width at half 
maximum (FWHM). Symbols are as in Fig. [3] The black solid 
line in the left panel indicates the observed correlation, whereas 
the black dashed line in the right panel shows the locus of equality. 
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Figure 5. As in Fig.|4]using for the line width the line dispersion. 



the intercept have not been significantly reduced with the 
enlarged sample. Similarly, the scatter stays unchanged at 
~ 50%, indicating that it is intrinsic rather than obser- 
vational. The scatter in the near-IR R-L relationship is 
slightly larger than that of the current optical R-L r elation - 
ship. For the entire optical sample of IBentz et alj (|2009a) 
of reverberation-mapped AGN with host galax y-subtracted 
fluxes (34 sources) we found a scatter of ~ 40% (|Landt et al.l 
l2011al . see their Table 2). 



5 THE NEAR-IR VIRIAL PRODUCT 

As we have shown in ILandt et all ^Ollbh . AGN black holes 
masses can be estimated from the near-IR virial product 
based on the 1 /im continuum luminosity and the width of 
the Pa/3 (or Paa) broad emission line because the widths 
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of the broad Paschen li nes are well-correla ted with those of 
the broad Balmer lines l|Landt et al,ll2008l ). The advantages 
of using the near-IR instead of the optical virial product 
are threefol d. First, since Paa and Pa/3 are observed to be 
unblended (|Landt et alj |2008). their width can be reliably 
measured. By contrast, the H/3 broad emission line is of- 
ten observed to be strongly blended with both Fe II and 
He II A4686 and often shows a "red shelf" most likely formed 
by weak Fe II multiplets and He I. Secondly, the AGN contin- 
uum around 1 fim is free from major contaminating compo- 
nents and can be easily determined, unlike the optical, which 
can suffer from a pseudo-continuum caused by blended and 
broadened Fe II emission. Finally, since the near-IR is much 
less affected by dust extinction, it can potentially be applied 
to dust-obscured AGN. 

In Figs. [4] and [5] (left panels) we show the relation be- 
tween the black hole mass derived from optical reverbera- 
tion mapping campaigns and the near-IR virial product. We 
have calculated the latter using for the line width the full 
width at half maximum (FWHM; Fig. |5| and the line dis- 
persion (aii nc ; Fig. [5| of the Paa or Pa/3 broad component, 
whichever line had the higher signal-to-noise ratio. Further- 
more, we have assumed that the R-L relationship has a log- 
arithmic slope of 0.5, since our results in Section[4]were con- 
sistent with this value. We have not plotted here the source 
Mrk 590, since it was found to have weak and noisy broad 
emission lines. A crucial step towards isolating the intrinsic 
broad-line profile is the subtraction of the narrow-line com- 
ponent. The narrow-line profile appears inflected, i.e., the 
transition point between the broad and narrow components 
is obvious, in roughly half of our sample (14/26 sources; 
see also Table Hand Fig. 0. Four sources (PG 0844+349, 
PG 0804+761, Ark 120, and NGC 3516) clearly lack a 
Paschen narrow-line component, since their profiles have a 
broad top. However, in the remainder (8/26 sources) the 
transition between the broad and narrow c o mpon ents is not 
perceptible. As discussed in lLandt et alj |2008), in these 
cases we have estimated the contribution of the narrow com- 
ponent to the total profile by fitting to its top part a Gaus- 
sian with FWHM equal to that of the narrow emission line 
[O III] A5007. In the case of the additional sample observed 
with GNIRS, for which we do not have contemporaneous 
optical spectroscopy, we have used instead the near-IR nar- 
row emission line [S III] A9531. This method assumes that 
the FWHM of [O III] A5007 or [S III] A9531 is representative 
of the narrow emission line region and subtracts the largest 
possible flux contribution from this region. 

As explained bv lPeterson et all (|2004T ). both measures 
of the line width have their advantages and disadvantages. 
Whereas the FWHM value can be measured directly in the 
spectrum and is less sensitive to line blending, the line dis- 
persion is well-defined for arbitrary line profiles and is less 
sensitive to the presence of even strong narrow-line com- 
ponents. In particular for the latter reason, the black hole 
masses from optical reverberation campaigns are always cal- 
culated using the line dispersion, whereby this quantity is 
measured in the rms spectrum, i.e., in the variable part of 
the spectrum. 

We have us ed the C routine MPFIT (version 1.1; 
lMarkwardtll2009l ). which solves the least-squares problem 
with the Levenberg-Marquardt technique, to fit both corre- 
lations. The best-fit lines (solid lines in Figs. [4] and [5]l are: 



log Mbh = (0.88 + 0.04) • (21ogFWHM + 0.51ogz/Li Mm ) 

-(17.39 + 1.02) (2) 

and 

log Mbh = (0.98 ± 0.04) ■ (2 log a nnc + 0.5 log vLi„ m ) 

-(20.02 + 1.11), (3) 

where Mbh is in solar masses, FWHM and o"ii ne are in 
km s _1 and vL\^ m is in erg s _1 . In both cases the corre- 
lation slopes are close to one, which suggests that the line 
width of an unblended broad emission line in a single-epoch 
spectrum is a suitable proxy for the line width measured in 
an rms spectrum. In order to verify this claim, we have plot- 
ted in Fig. [7] the line width of the H/3 broad emission line as 
measured in the rms spectrum versus that of the Paa and 
Pa/3 broad components for both cases, i.e., using the FWHM 
(upper panel) and the line dispersion (lower panel). In the 
case of FWHM, we find a correlation slope consistent with 
unity of 1.2 ± 0.1, however, the line dispersion in rms spec- 
tra appears to be smaller than that in single-epoch near-IR 
spectra (correlation slope of 0.5 + 0.1). Since the line disper- 
sion strongly depends on the total flux of the emission line, 
this result could be explained if the broad wings vary less 
than the core of the emission line since they are composed 
of optically thin gas (jShields et al.|[l995l '). 

In order to estimate black hole masses from single-epoch 
near-IR spectra one then has two options: 1. directly apply- 
ing the relationship between black hole mass and near-IR 
virial product (eqs. ([2)| and (0); and 2. applying the near- 
IR R-L relationship together with the virial theorem (eq. 
HJ). In Figs. |3]and[S] (right panels) we show the results for 
the latter option using scaling factor values of / = 1.4 and 
5.5 in the case of FWHM and line dispersion, respectively 
llOnken et al 1I2004D . Clearly, both measures of the line width 
give results consistent with the reverberation-based black 
hole masses. But we note that in the case of the line dis- 
persion the second (but not the first) option is expected to 
overestimate the black hole mass (by a factor of < 2) given 
the relationship between the H/3 rms line width and that of 
the Paa and Pa/3 broad components presented in Fig. [7] 



6 LOW-LUMINOSITY AGN 

Four AGN in our sample, namely, Mrk 590, NGC 3227, 
NGC 3516 and NGC 4051, were found to be in a very low 
luminosity state with a continuum SED domin ated by host 
galaxy emission (see Fig. [2] top panels, and lLandt et alj 
(|2011br i. their Fig. 6). Therefore, in these cases the ob- 
served 1 /jm continuum luminosity represents only an up- 
per limit to the AGN value. As Fig. [3] shows, these sources 
(shown as green filled circles) populate indeed the low end of 
the near-IR luminosity distribution of our sample. However, 
they follow not only the near-IR R-L relationship remark- 
ably well, but also t he so-called black hole mass-luminosit y 
(M-L) relationship (|Kaspi et al.ll2000l : |Peterson et al.|[20o3 l ■ 
which we plot in Fig. [8] for the near-IR. Therefore, the AGN 
1 /im continuum luminosity of these sources cannot be sig- 
nificantly (more than a factor of ~ 2 — 3) lower than that of 
the host galaxy. 
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Figure 6. Pact or Pa/3 broad-emission line profiles in velocity space relative to the expected rest frame wavelength (thick lines). The 
profiles have been continuum-subtracted and normalized to the same peak intensity (of the broad component). The subtracted narrow 
components are also shown (thin lines). 



Can and should black hole masses be estimated for AGN 
with 1 fim continuum luminosities dominated by host-galaxy 
starlight? The near-IR virial products and near-IR black 
hole mass values calculated for the four host-galaxy domi- 
nated sources in our sample, although only upper limits, fol- 
low well the general relationships (see Figs. [4] and [5|. But we 
also observe that the SEDs of all these sources are dominated 
by the hot dust component at wavelengths A J> 1.5 fim, i.e., 
in the near-IR H and K bands. 



7 SUMMARY AND CONCLUSIONS 

We have presented an improved calibration of our near-IR 
relationship for estimating black hole masses by increasing 
the sample of reverberation-mapped AGN with high-quality 
near-IR spectroscopy from 14 to 23 sources. The additional 
observations were obtained at the Gemini North observatory 
with the Gemini Near-Infrared Spectrograph (GNIRS). Our 
main results can be summarized as follows. 

(i) The near-IR radius - luminosity (R-L) relationship 
has in most cases a best-fitting slope of ~ 0.6 ± 0.1. This 
value is slightly steeper than our previous result but consis- 
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Figure 7. The width of the H/3 broad emission line as measured in 
the rms spectrum versus that of the Paa or Pa/3 broad component, 
whichever line had the higher signal-to-noise ratio. The line width 
is measured as the full width at half maximum (FWHM; upper 
panel) and the line dispersion (lower panel). Symbols are as in 
Fig. [3] The solid lines show the observed correlations and the 
dashed lines mark the locus of equality. 



tent with a value of 0.5, which is expected based on simple 
photoionization arguments. The scatter in the relationship 
stays unchanged at a value of ~ 50%, indicating that it is 
intrinsic rather than observational. The additional sample 
has improved the number statistics in particular at the high 
luminosity end. 

(ii) The black hole mass derived from optical reverber- 
ation mapping campaigns correlates strongly with the near- 
IR virial product, i.e., the product between the integrated 
total 1 fira continuum luminosity and the width of the un- 
blended Paschen hydrogen broad emission lines Paa and 
Pa/3. The correlation slope is ~ 1, which suggests that the 
line width of an unblended broad emission line in a single- 
epoch spectrum is a suitable proxy for the line width mea- 
sured in an rms spectrum. 

(iii) We have excluded four AGN in our sample from 
the calibration of the near-IR relationship, since they were 
found to be in a very low luminosity state and so to have 
a continuum emission dominated by host galaxy starlight. 
Nevertheless, these sources still follow the near-IR relation- 
ship remarkably well, indicating that it can be used to re- 
liably derive upper limits on the black hole masses of low- 
luminosity AGN. 

In the future we plan to start a near-IR reverberation 
mapping campaign, which will allow us to obtain time lags 
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Figure 8. Black hole mass determined from optical reverberation 
campaigns versus the continuum luminosity at 1 fim. Symbols are 
as in Fig. [3] 



for the Paschen broad emission lines and thereby black hole 
mass relationships fully calibrated in the near-IR. 
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APPENDIX A: GEMINI NEAR-IR SPECTRA 

This paper has been typeset from a TjtX/ DTjtX file prepared 
by the author. 
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PG 0052 + 251 (z = 0.155) 
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Figu re Al. Gem ini GNIRS near-IR spectra shown as observed flux versus rest-frame wavelength. Emission lines listed in Table 4 of 
lLandt et alj ||200S|) are marked by dotted lines and labeled; black: permitted transitions, green: permitted Fe II multiplets (not labeled), 
red: forbidden transitions and cyan: forbidden transitions of iron (those of [Fe II] not labeled). 
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PG 0026+129 (z = 0.145) 
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PG 0804 + 761 (z = 0.100) 
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Figure Al - continued 



© 2013 RAS, MNRAS 000.ITHTT1 



H. Landt et al. 





I 






1 1 
























1 ■ 1 










a 
















CM 










S "S 


I 














X 


EE 




o 




o 


"E 


Y 


Pfi 1 


Pn 1 


(£ 




a 


;; cl 






o 














5 






















I 1 




Hel 














, 1 






M 1 












r 












1 











PG 1211 + 143 (z = 0.081) 
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Figure Al - continued 
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Figure Al - continued 
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NGC 3516 (z = 0.009) 
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Figure Al - continued 
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